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ABSTRACT 

We propose to study the behavior of complicated numerical solutions to Einstein's 
equations for generic cosmologies by following the geodesic motion of a swarm of test 
particles. As an example, we consider a cylinder of test particles initially at rest in the 
plane symmetric Gowdy universe on T 3 x R. For a circle of test particles in the symmetry 
plane, the geodesic equations predict evolution of the circle into distortions and rotations 
of an ellipse as well as motion perpendicular to the plane. The evolutionary sequence of 
ellipses depends on the initial position of the circle of particles. We display snapshots of 
the evolution of the cylinder. 



I. Introduction. 

Modern computer technology allows numerical investigation of nontrivial solutions to 
Einstein's equations. In particular, a program of research is now in progress to discover 
the nature of the Big Bang in generic spatially inhomogeneous universes [1]. Yet this 
ability to tackle complicated spacetimes brings with it the problem of how to interpret the 
computer's output. This problem is especially acute in the cosmo logical context where the 
interpretational familiarity of an asymptotically flat region is absent. 

As an example, consider spatially inhomogeneous vacuum cosmologies on T 3 x R that 
possess a single spatial U(l) symmetry (the U(l) models). According to Moncrief [2, 1], 
these models can be described by a metric written in terms of five degrees of freedom at 
each spatial point. (These are equivalent to the two degrees of freedom that describe the 
gravitational waves with the others restricted by constraints and coordinate conditions.) 
Clearly, graphical display of the amplitudes of all the variables on the spacetime coordinate 
grid will not yield a physical understanding of the model. We propose as an alternative a 
numerical gravitational wave detector: the geodesic motion of a swarm of test masses in 
the computer generated spacetime. 

As an initial implementation of this idea, we present here snapshots from computer 
movies of the geodesic motion of a cylinder of test masses in the simpler Gowdy T 3 x R 
cosmology [3]. (Although we shall continue to refer to a cylinder, T 3 topology means 
that the detector is actually a 2-torus.) The numerical methodology for solving Einstein's 
equations for this model (described elsewhere [1]) allows it to be regarded as a test case for 
the U(l) models. The Gowdy model can be interpreted as plane symmetric gravitational 
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waves propagating in one direction in a background cosmology with the two polarizations 
of the waves nonlinearly coupled. Extensive numerical investigations [1, 4] have shown 
that as Gowdy models evolve toward the singularity, the nonlinear wave interactions cause 
the appearance of ever shorter wavelength modes. Depending on initial conditions, this 
nonlinear growth gets suppressed as the model becomes velocity dominated [1, 5]. (This 
means that time derivatives dominate spatial derivatives in Einstein's equations so that 
each point in the universe evolves independently causing the spatial profile of the metric 
variables to become frozen.) 

For a particular class of initial conditions, we shall compare the behavior of our "de- 
tector" in a universe with little growth of small scale spatial structure to that in one with 
significant growth. We shall consider both "swarm of particles" and "rendered spatial 
volume" visualizations of the detector. These, in turn, will be compared to the amplitudes 
of the metric variables. Directions for future enhancements of the detector will also be 
discussed. 

II. The Gowdy T 3 x R Model. 

The metric for this cosmology is 

ds 2 = e 2a {-dt 2 + dO 2 ) + te p [da 2 + (e~ 2P + Q 2 ) d5 2 + 2Q da dS] (1) 

where P, Q, and a depend only on 9 and t. The two Killing fields are (d /da) and (d/dd). 
Einstein's equations for this model consist of wave equations for P and Q: 

P, TT - e- 2T P, 0e - e 2P (Q 2 - e~ 2T Q 2 e ) = 0, (2a) 
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<2,rr - e- 2T Q m + 2 (P, T Q, T - e~ 2T P, Q, e ) = 0, 



(26) 



and the momentum and Hamiltonian constraints respectively for A in terms of P and Q: 



X, e -2(P, e P, T +e 2P Q,eQ,r)=0, 



(3a) 



A, T - K + e-^P 2 + e 2P (Q 2 + Q* )] = 



(36) 



where X = 4a + t and r = — Int. 
III. Geodesic equations. 

In order to construct our gravitational wave detector, we require the geodesic equations 
for the Gowdy spacetime. It is easily seen that symmetries allow the detector to have a 
simple shape. For u a the tangent vector to a geodesic and £ a a Killing vector, u a ^ a is a 
constant of the motion. It then follows that the geodesies in the Gowdy spacetime admit 
the following two constants of the motion: 



ci = t e 



da ^ dS 
— + Q — 
dp dp 



(4a) 



c 2 = t e 



~ da / -op d>5 
Q + (e + g.) _ 



(46) 



where p is the affine parameter of the geodesic. In particular, geodesies for which da / dp 
and dS/dp vanish at the initial time, have c\ = c-i = so that a and 8 are constant. The 
motion of the particle is thus given by finding 6 as a function of t. This means that the 
detector can be constructed as a cylinder of test masses with fixed coordinates in the a-S 
plane. 



Note that, for the geodesies that describe the detector, the proper time AT is given 



by 



/ 



AT = / Ldt 



(5) 



where the Lagrangian L is given by 



L = e a 1 



d6 



1/2 



(6) 



Since geodesies are those paths that extremize proper time, the geodesic equation reduces 
to the Euler-Lagrange equations for the Lagrangian L, 



f9_ 

dt 2 



+ 1 



de 

~dl 



' ' / da dO da \ 

\dt dt + de) = 



(7) 



Note that, generically, each ring composing the cylindrical detector must move in 9 as the 
spacetime evolves. 

In studying the approach to the singularity, it is helpful to use the coordinate r since 
the singularity occurs at r = oo. In terms of this time coordinate the equation of motion 
becomes 

rftO rift ( . T^l 2 \ / _ Fin rift f)n\ 

(8) 



dr 2 



d6 

Tr + 



,2t 



dr 



_ 2r da dO da\ 
dr dr 
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Our approach then is to integrate equation (7) or (8) for representative rings of the cylin- 
drical detector in the numerically generated background spacetime (see Fig. 1). 
IV. Motion of a ring of particles. 

Consider one of these rings with the given initial values of (0,t,dO/dp,dt/dp). We 
wish to follow the shape of the ring as the particles move. The motion of the ring of test 
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particles is described by a single solution of equation (7); the particles which compose 
the ring differ from each other only in the (constant) values of the coordinates a and 5. 
However, as the particles move they encounter changing metric functions P and Q. This 
gives rise to changing lengths associated with the constant spatial coordinates a and 5 
and thus a changing shape of the ring. From equation (1), the conformal metric of the 
symmetry plane is given by 

d s 2 = e p d(j 2 + (e~ p + e p Q 2 ) d5 2 + 2e p Qdad5 . (9) 

Since all particles in the ring are described by the same 6{t), distances between the particles 
are given by \ftds. The overall size of the ring is then proportional to \ft. However, since 
we are only concerned with the shape of the ring, we will consider only distances measured 
by the conformal metric ds 2 . 

We first treat the simpler case of a polarized Gowdy spacetime (Q = 0) where 

d~s 2 = e p da 2 + e~ p d5 2 . (10) 

In terms of X = e p / 2 a and Y = e~ p l 2 5 : the conformal metric is flat. Suppose that at the 
initial time the ring has the shape of a unit circle: Xq + Y 2 = 1. This means that 

e Po a 2 + e -Po 6 2 = 1 (n) 

Since the (a, 5) coordinates of the particles are constant, equation (11) holds at all times. 
Expressing (cr, 8) in terms of (X, Y) , we find 

e Po~P X 2 + e P-P« Y 2 = l ^ (12) 
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the equation for an ellipse. Thus the ring of particles that started out as a circle is distorted 
into an ellipse by the gravitational waves. 

We now consider the general case of an unpolarized Gowdy spacetime. Define the 
coordinates X and Y by X = e p / 2 (a + QS), Y = e~ p / 2 5 yielding ds 2 = dX 2 + dY 2 . 
However, the X axis rotates relative to a parallel propagated frame. Let ip be the angle of 
rotation and define (X, Y) by 

X = X costfj + Y sin^ , (13a) 

Y = Y cos^ - X simf; . (136) 

Then ds 2 = dX 2 + dY 2 and the X axis is nonrotating. Thus expressing the positions of 
the particles in the ring in terms of the (X, Y) coordinate system will tell us how much 
the ring is rotated as well as how much it is distorted. 

We now calculate the rotation angle ip. Let (n a ,m a ) be an orthonormal basis in the 
(X,Y) directions (described e.g. in the nonrotating coordinate system). Let u a be the 
tangent vector to the geodesic at the center of the ring. Then the changes in the basis 
vectors are related by 

u a V a n b = (u a V a ^) m b . (14) 
For £ a = (d/da) a , the metric and the definition of (X, Y) lead to 



i a = te p (d a a + Qd a 5) 
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(16) 



Using the fact that £ a is a Killing vector, we find 

u a V a n b = 6 u a d a (V 1 / 2 e" p / 2 ) + e~ p / 2 u a d [a ^ h] . (17) 

Some straightforward algebra yields 

u a V a n h =^ 2 e p ' 2 (u a V a Q) d b 5 
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We therefore have 



1 e p (u a V a Q) m h . (18) 



u a V a i(j = ^e p u a V a Q 



and thus 



1 f P fdQ dQ d9\ 



We can now find the shape of the ring of particles. Since the ring starts out as a unit 
circle, we have 

1 = X 2 + Y 2 = e p » (a + Q 5f + e" p ° 5 2 . (20) 

Since a and 5 are constant, equation (20) holds at all times. Expressing a and 5 in terms 
of X and Y gives 

e p °~ p [X + (Q - Q) e p Y] 2 + e p ~ p "Y 2 = 1 . (21) 



This is again the equation of an ellipse. Finding the angle ip then allows us to rotate this 
ellipse to the (X, Y) coordinate system. 
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Our explicit procedure is to construct an initial set of unit circle rings uniformly 
distributed along the 6 axis with zero initial velocity in 9. A previously computed numerical 
spacetime {P, Q, a}(#, r) is used in the geodesic equation (8) to evolve 6>(r) for each ring. 
When preselected values of r are reached, the shape of each ring is computed using equation 
(21) with equations (13) to construct the ellipse in a nonrotating frame. For example, 
expressing (X, Y) in terms of polar coordinates (r, () in the a-5 plane gives r(() for the 
ellipse directly from equation (21). 

V. Results. 

Various techniques have been tested to visualize the evolution of the detector. Two 
of these are illustrated in the figures. First we consider explicit display of the detector's 
component rings. (The visualization here was constructed using D 2 Software MacSpin 
3.0 on a Macintosh Quadra 660AV. A similar construction using custom software written 
by one of us (VS) runs under X-Windows on a Unix workstation.) Alternatively, we 
can display a rendered volume of the entire cylinder with interpolations to find values in 
between the selected rings. The displayed values are defined as the region inside each ring. 
The location of the ellipse in a 4 x 4 square was found by transforming each point of an 
array on the square to the coordinate system defined by the major and minor axes of the 
ellipse. (This visualization uses Spyglass Dicer 3.0 on the Macintosh. One of us (VS) has 
written custom software for visualization under X- Windows.) 

The detector at r = is shown in Fig. 1 as a unit cylinder. Figure 2 illustrates the 
distortion of the ellipses for a polarized Gowdy background later in the evolution. Note 
that, as expected, there is no rotation of the ellipse. Agreement with equation (11) has 
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been checked. 

More interesting is a comparison between two Gowdy universes which differ only in the 
value of a parameter v$ in their initial conditions. This parameter controls the nonlinear 
growth of small scale spatial structure in these models [1, 4]. In Figs. 3a and 3b, the 
complete evolution of P, Q, and A is displayed for vq = 0.5 and vq = h respectively. 
Note the increased spatial structure in P for the larger v . (The interpretation of these 
features is discussed elsewhere [1, 4].) Figures 4 and 5 display snapshots of the evolution 
of the detector (using both visualizations) at three values of r for v$ = 0.5 and vo = 5 
respectively. In all cases, the values of the metric functions P and Q (which determine the 
detector's shape in the a-5 plane) are also displayed. Within the limited resolution of the 
detector as illustrated here (due to convenience rather than principle), the features of the 
spacetime are clearly seen. Note both the spatial variation in orientations of the ellipses 
and the rotation (as well as distortions) of the ellipses with time. 

V. Discussion. 

We have demonstrated the usefulness of a numerical gravitational wave detector in the 
visualization of spatially inhomogeneous cosmologies. For application to the U(l) models, 
the absence of the symmetry plane probably means that the detector should consist of a 
swarm of individually followed test particles. The initial conditions of the swarm could 
be adjusted to explore a variety of features of the U(l) spacetime on a variety of spatial 
scales. 
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FIGURE CAPTIONS 



Fig. 1. The initial state of the detector, (a) The volume of space within the unit rings in 
the a-5 plane that compose the detector. The a-5 plane is represented as a 25 x 25 array 
of points for — 2 < X, Y < 2. (See Section IV.) (b) Each of 10 rings that constitutes the 
cylinder is displayed as 100 points at fixed angle Q in the a-5 plane. In each case, 10 frames 
of the evolution are saved. 

Fig. 2. The detector in a polarized model. A polarized model with P(6>, 0) = 0, P, T (6>, 0) = 
0.5 cos 6>, Q(6>, 0) = 0, and Q, T (6>, 0) = is evolved up to r = 9.2. The shape of the detector 
at r = 5.2 is displayed for both visualizations. 

Fig. 3. The unpolarized background spacetimes. Unpolarized models with P(6>, 0) = 0, 
P, T (9, 0) = v cos#, Q(9, 0) = cos#, and Q, r (9, 0) = are evolved up to T f = 8.40. The 
complete evolution of the metric functions P, Q, and A is shown for (a) vq = 0.5 and (b) 
v = 5. In all cases, < 9 < 2n and < r < 77. The vertical scale indicates the value of 
the appropriate metric variable. 

Fig. 4. Snapshots of the detector's evolution in the unpolarized model with vq = 0.5. 
Both visualizations of the detector and slices in time of P and Q vs. 9 are shown for (a) 
r = 2.80, (b) r = 4.67, and (c) r = 6.53. For the metric variables, the horizontal scale is 
< 9 < 2n while the vertical scale is the same as that in Fig. 3a. 

Fig. 5. The same as Fig. 4 for vq = 5. The vertical scale for the metric variables is from 
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Fig. 3b. Here 20 rings were used to describe the detector. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9405069vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9405069vl 



This figure "fig3-l.png" is available in "png" format from: 
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